Did Boomerang hit MOND ? 
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Purely baryonic dark matter dominated models like MOND based on modification of Newtonian 
gravity have been successfull in reproducing some dynamical properties of galaxies. More recently, 
a relativistic formulation of MOND proposed by Bekenstein seems to agree with cosmological large 
scale structure formation. In this work, we revise the agreement of MOND with observations in light 
of the new results on the Cosmic Microwave Anisotropies provided by the 2003 flight of Boomerang. 
The measurements of the height of the third acoustic peak, provided by several small scale CMB 
experiments have reached enough sensitivity to severely constrain models without cold dark matter. 
Assuming that acoustic peak structure in the CMB is unchanged and that local measurements of 
the Hubble constant can be applied, we find that the cold dark matter is strongly favoured with 
Bayesian probability ratio of about one in two hundred. 



I. INTRODUCTION 

The measurements of the Cosmic Microwave Back- 
ground (CMB) anisotropies, most notably by the 
Wilkinson Microwave Anisotropy Probe (WMAP) 
mission have truly marked the beginning of the 
era of precision cosmology. In particular, the shape 
of the measured temperature and polarisation angu- 
lar power spectra are in spectacular agreement with 
the expectations of the standard model of structure 
formation, based on primordial adiabatic and nearly 
scale invariant perturbations (see e.g. 0). More re- 
cently, ground based and balloon borne CMB experi- 
ments like extended VSA 0LCBI [|, Acbar [|, DASI 
@ and BOOMERANG-03 0, have probed the CMB 
power spectra at smaller scales, confirming the pres- 
ence of acoustic oscillations and providing the first 
unambiguous detection of polarisation. 

Moreover, new, complementary, results from the 
Sloan Digital Sky Survey (SDSS) on galaxy clustering 
(see e.g. ||) and, more recently, on Lyman-ev Forest 
clouds [!j are now further constraining the scenario. 

Since all these measurements appear in spectacular 
agreement with the ACDM model, based on a cosmo- 
logical constant and on cold dark matter, is definitely 
timely to investigate what space is left for alternative 
theories. 

Perhaps the most exotic alternative to the standard 
model one could consider is MOdified Newtonian Dy- 
namics (MOND, 0]) where a purely baryonic model 
with modifications to standard (Newtonian) gravity is 
suggested. In MOND the departure from Newtonian 
law a = — V<I>n is given by: 

A(|a|/o )a= -V$ N (1) 

where $n is the Newtonian potential of the visible 



matter, do is an acceleration scale while fl(x) ~ x for 
x < 1 and ji(x) -> 1 for x > 1. If o « 1 x 1CT 8 
cm s -2 the Newtonian law is recovered in the solar 
system where accelerations are large compared to do- 
The above empirical formula has been originally 
proposed to explain the fact that rotation curves of 
disk galaxies become flat outside their central parts. 
While in the standard dark matter paradigm flat rota- 
tion curves are explained by assuming a spherical halo 
of invisible dark matter around visible disk galaxies, in 
MOND there is no need to include non baryonic dark 
matter since, thanks to the Eq^ galaxies far out ex- 
hibit an approximately spherical Newtonian potential 
without the inclusion of the dark matter. Attempts 
were made to confront MOND with clusters of galax- 
ies (see e.g. 0,0) and the large scale structure (e.g. 
[l3|) with mixed success. 

While the non-baryonic dark matter paradigm is 
definitely more compelling for its aesthetic simplicity 
than a modification to Newtonian gravity, the MOND 
model has been claimed successful in other aspects 
(see e.g. 0,^|) and MOND proponents insist that 
this alternative model for gravity merits serious ex- 
amination. 

The MOND theory has suffered from a lack of a 
successful relativistic formulation, that would allow 
one to compare it to observations of CMB and Large 
Scale Structure. Nevertheless, there were some at- 
tempts to confront CMB data 16], which find that 
the first 2 observed acoustic peaks in the CMB spec- 
trum are compatible with MOND at the price of a 
substantial neutrino mass, which is barely compatible 
with current laborator y b ounds |l7j| . or at the price 
of including curvature [18| , which is at odds with the 
inflationary scenario. 

A major step in the direction of developing the sim- 
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pie MOND formula into a more robust theory of grav- 
ity has been recently proposed by Bekenstein [Igj . 
In this paper, a relativistic gravitational theory has 
been presented whose nonrelativistic weak accelera- 
tion limit accords with MOND while its nonrelativis- 
tic strong acceleration regime is Newtonian. 

Moreover, Bekenstein's model provides a specific 
formalism for constructing cosmological models and 
testing MOND using cosmological data. Indeed, more 
recently, Skordis et al. [20j produced the first theoret- 
ical prediction for CMB anisotropies and Large Scale 
Structure in the case of Bekenstein's model. It has 
been shown that the Bekenstein model may be put 
in agreement with the WMAP data and Large Scale 
Structure observations. Similar to previous results au- 
thors find agreement if neutrinos ensure that peak po- 
sitions are unchanged. The results are obviously of 
great relevance since the model has no cold dark mat- 
ter in it and may therefore be considered as an impor- 
tant alternative to the present cosmological scenario, 
which assumes a fine-tuned cosmological constant and 
yet to be discovered dark matter particles. 

In this brief report we point out that recent, small 
scale, CMB data already provide discriminating power 
between these two scenarios. In the standard cold 
dark matter scenario, the amplitude of the CMB peaks 
is sensitive to the amount of dark matter because of 
two effects: increasing the matter density on one hand 
decreases the radiation driving while on the other 
hand it increases the depth of potential wells. These 
two effects nearly cancel out in the amplitude of the 
second acoustic peak, but conspire to produce a higher 
amplitude of the third peak (see |2l|)- The height of 
the second peak to the first peak therefore contains 
information on the baryonic content of the Universe, 
while the ratio of the third peak to the first peak 
height tells us about the matter density. 

A generic prediction of purely baryonic dominated 
models like MOND is therefore that peaks in the CMB 
power spectrum should be strictly decreasing in am- 
plitude. 

Extraordinary WMAP results on the first two 
peaks, coupled with the recent small scale CMB data 
on the third peak now have enough power to discrim- 
inate between these two scenarios and to determine 
the amount of cold dark matter. 

The goal of this brief report is to examine the data 
on the third peak, with special emphasis on the re- 
cent measurements of the CMB fluctuations by the 
Boomerang experiment. 

The ultimate test of MOND, would be a full con- 
frontation of the relativistic theory with the data. 
This is a daunting task, given that the theory is very 
complicated with yet to be fully understood pertur- 
bation theory and several free parameters including 
a free function. We note, however, that models dis- 
cussed in |2(j are in a complete agreement with generic 
A-CDM predictions in the range (. > 200. We there- 
fore assume that this is a generic prediction of the 
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Prior 



Vfr: 



log lO 10 ^ 



^run 



BBN (see text) 
(0,0.99) 
(0.5,10) 

(0,1) 
(0.5,1.5) 

(2,5) 
(-0.2,0.2) 



TABLE I: Flat priors on the cosmological parameters. No- 
tation (x,y) implies a a flat prior between x and y. Run- 
ning spectral index (last parameter) was used only in the 
part discussed in the last paragraph of section ITTT1 



Bekenstein's theory and proceed by fitting the £ > 200 
region of the CMB data with the standard A-CDM 
models to see whether models with zero cold dark mat- 
ter density are compatible with the data. Whether our 
assumption is a justified one is to be seen, however, we 
feel it nevertheless provides a first order confrontation 
of the data with the theory. Our approach is orthog- 
onal to that of [2(j in the sense that it provides con- 
straints on any theory that leaves the CMB physics 
unchanged on scales smaller or roughly equal to that 
of the first acoustic peak. 



II. ANALYSIS 



We use the Cosmo-MC package ( [23] ) to perform pa- 
rameter estimation on standard flat A-CDM models 
using top- hat priors on the following 6 parameters: 
uj\y = Q,^h 2 (the baryon density of the universe), Wdm 
(the dark matter density of the universe), 9 (ratio of 
the sound horizon to the angular diameter distance 
to the surface of last scattering, multiplied by 100), 
^frac (the fraction of dark matter in form of massive 
neutrinos), n s (spectral index of primordial fluctua- 
tions), log A s (the amplitude of primordial scalar fluc- 
tuations). The priors used are listed in the Table [I] 

In our parametrisation the density of the cold dark 
matter is given by 



^cdm — 



h 2 



(2) 



We intentionally omitted r, the optical depth to 
the last scattering, from our parametrisation as it is 
completely degenerate with the amplitude in the mul- 
tipoles of interest. 

The following datasets were used in our parametri- 
sation: WMAP 0, E|, VSA 0], CBI H, Acbar 
0, and the latest Boomerang results Q. In these 
datasets we have used the default Cosmo-MC distribu- 
tion datasets for VSA, CBI and Acbar experiments. 
In all datasets any points with I < 200 were removed 
and additionally all points with £ < 375 were removed 
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FIG. 2: This plot shows the marginalised probabil- 
ity distribution for f2 C md for the WMAP dataset alone 
(thin line) and WMAP / VS A / ACB AR/CBI /Boomerang 
datasets (thick line). Curves are normalised to the same 
total area under the curve. 



FIG. 1: This plot shows the experimental data used and a 
few theoretical spectra. Points with large error-bars were 
removed from the plot (but actually used in the param- 
eter estimation chains) to maintain clarity. The exper- 
iments are as follows: WMAP (bars), VSA (triangles), 
CBI (crosses), ACBAR (crossed circles), Boomerang (cir- 
cles) . The theoretical graphs plotted are most likely model 
in our chains for the all CMB data (solid line) and two 
models picked from best fit models with Q m < 0.01 from 
each chain sampling WMAP data only. Parameters for 
dashed model are (wi,,cj v , h, n s ) = (0.22,0.096,0.76,0.59) 
and those for the dotted are (0.22,0.15,0.78,0.43). See 
text for discussion. 



from Boomerang temperature data (to prevent cosmic 
variance coupling to the WMAP dataset). 

The basis of our analysis are the heavily over- 
sampled chains of the WMAP dataset (with about 
^260,000 independent points) on the parametrisation 
described above, together with a Big Bang Nucleosyn- 
thesis (BBN) prior of uo b = 0.020 ±0.002 Q. There is 
an additional prior on 0.4 < ft, < 1.0 hard-coded into 
the Cosmo-MC code and our posterior space is cut by 
the lower end of this prior. Releasing this prior would 
weaken our constraints, but h < 0.4 models would 
face many difficulties with other cosmological probes. 

We importance sample these chains with additional 
data in order to get improved constraints. Apart from 
savings in the CPU time, the importance sampled 
chains still include models from the low Od m region 
which is of interest for the MOND models and would 
otherwise be present only far into the tails of the prob- 
ability distribution. We also made a similar run with- 
out the BBN prior but with very wide flat prior on Wb 
instead, to check the effect of BBN prior dependence. 



III. RESULTS 

The parameter of major interest here is Q c dm- A 
MOND predicts a zero f2 c dm, while the standard pre- 
diction is that of a O c( jm of about 0.3. In the Figure Q] 



we plot the data we used (note that some points were 
omitted from the plot - see caption) and a few theo- 
retical predictions. The solid line model corresponds 
to a standard ACDM flat model and fits the data very 
well. The dashed and dotted models were cherry- 
picked from the models that have the highest likeli- 
hood of models that satisfy Vl m < 0.01 in each indi- 
vidual MCMC chain (using WMAP data alone). The 
dashed model illustrates the common wisdom that it 
is possible to construct models that have nearly iden- 
tical peak positions and heights of even peaks with 
zero dark matter. Finally, the dotted model shows 
an example of a zero dark matter model that is al- 
lowed by the WMAP data but obviously at odds with 
measurements of the small scale power. 

In the Figure [21 we plot the marginalised probabil- 
ity distribution for fi c d m for the WMAP dataset and 
the WMAP dataset after inclusion of all the other 
CMB data. We note that the WMAP alone admits 
the zero cold dark matter solutions, in agreement with 
previous investigations 0, |2(J . The addition of other 
datasets, however, strongly rejects this region of the 
parameter space, without resorting to non-CMB ex- 
periments. 

In order to quantify this further, we use two sta- 
tistical tests. Firstly, we compare two basic models, 
namely the flat CDM model with £! c d m between and 
1 and the MOND model for which f2 c dm = using 
Bayesian model comparison pEj . In order to do this 
we compute the logarithm of evidence ratio with a 
version of the Savage-Dickey test, equivalent to that 
described in the Appendix of (2(|. Secondly we esti- 
mate the number of "sigmas" in a frcquentist manner 
by comparing the likelihoods of the f2 c dm = point 
and the most likely point in the marginalised probabil- 
ity distribution for f2 c dm- The likelihood ratio can be 
converted to a number of standard deviations using 

the prescription n n = < / 2 „ £max , which returns the 
* 1 V -Cmond ' 

expected result in the Caussian case. In such short 
communication it is impossible to compare the two 
methods in depth, but we note that frequenstist ap- 
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Dataset 


A logS 


n a 


WMAP 


0.1 


1.5 


WMAP / VS A / ACB AR/CBI 


2.1 


2.6 


WMAP/VSA/ACBAR/CBI/BOOM 


3.6 


3.1 


4-HST 


5.2 


3.6 


+HST +SN 


10.5 


5.1 


ALL CMB w/o BBN 


2.2 


2.6 


-BBN +HST 


3.2 


3.3 


-BBN +HST +SN 


4.3 


3.6 



TABLE II: This table shows results of the two statistical 
tests described in the text for a range of datasets con- 
sidered. The sign convention is such that higher number 
implies higher statistical evidence in favour of cold dark 
matter models. The AlogE is trivially interpreted as the 
logarithm of probability ratio. See text for discussion. 

preaches neary always give higher "detection" confi- 
dences. Bayesian evidence ratio has an advantage that 
it directly encodes the probabiltiy ratio. 

The results are summarised in the Table|H] The ex- 
act numbers somewhat depend on the binning width 
and therefore the numbers in the table are accurate to 
about 0.1 in both columns. Some row state the con- 
straints upon adding two extra constraints on top of 
all CMB data. HST label denotes the Hubble Space 
Telescope (HST) constraint on the Hubble's constant 
|27j . which is conveniently described as a Gaussian 
around h = iJ /100km/s/Mpc = .72 with 1-er dis- 
persion of 0.8. The HST data actually has a rather 
strong effect on our results as the Hubble constant 
favoured by the fi c( j m = models is rather low. SN 
label denotes the gold dataset of the |23| supernovae 
data. These results must be taken with some cau- 
tion, because it is not entirely clear that the standard 
interpretation of these two cosmological probes is ap- 
plicable in the MONDian setting. 

Finally, we have also considered running of the spec- 
tral index, defined by n run = dn s /dlnfc with pivot 
scale set to k = 0.05 Mpc/i -1 . Using a top-hat prior 
between -0.2 and 0.2 on n run we find that this param- 
eter is completely unconstrained by the WMAP data 
and only weakly constrained by the all CMB data. 
However, the very high I points from Acbar seem to 
favour negative running indices, thus even lowering 
the third peak required by MOND. Consequently we 
find that the bound on fl c dm is unaffected when all 
CMB data are included. 



IV. DISCUSSION AND CONCLUSIONS 

In this brief report we have analysed the latest CMB 
data in light of the recently renewed interest in the 
MOND theories of gravity. We simplified (and poten- 
tially oversimplified) the theory by assuming that the 
small scale CMB fluctuations are unmodified by the 



MOND theory in accordance with recent attempt [2(i| 
to model linear fluctuations in the relativistic theory 
of MOND recently proposed by Bckcnstcin. Under 
this assumption we find that the WMAP data alone is 
fully consistent with the fldm = required by MOND, 
in accordance with previous findings. A component 
of massive sterile neutrinos is required with ui v > 0.1, 
which is marginally consistent with earth-based beta 
electron decay experiments. Addition of other CMB 
data constrains the third peak height, which encodes 
the information on the presence of the cold dark mat- 
ter. The data before the latest Boomerang dataset 
weakly favour the cold dark matter but only at 1 in 8 
probability ratio for a Bayesian analysis and less than 
3-<7 for the more conventional approach. The VSA's 
measurement of power at the top of the third peak is 
probably crucial in absence of Boomerang data. Ad- 
dition of the third peak from latest Boomerang data 
gives crucial extra information. The Bayesian proba- 
bility ratio in favour of cold dark matter models in- 
creases to 1 in 36, while the likelihood ratio breaks the 
3-<7 "barrier" . We note that using less general models 
for the CDM setting would make the Bayesian model 
selection even stronger in favour of cold dark matter 
models. 

Releasing the BBN prior somewhat weakens the 
constraints. We note however, that there is no ob- 
vious mechanism how could a MOND theory evade 
BBN constraints and that such models are addition- 
ally characterised by very low values of h and n s . 
Marginalised value of to h ~ 0.025 for these models. 

Finally, we have added two other standard cosmo- 
logical probes, the HST measurements of the Hubble's 
constant and the recent supernovae measurements of 
the luminosity distance. Taken at the face value, they 
seem to blow the MOND model into oblivion. Care 
must be taken, however, in interpreting these two 
datasets as it is not clear whether it is appropriate to 
include them in the MONDian scenarios without de- 
tailed treatment of possible MONDian effects on the 
background evolution. 

We have shown that running of the spectral index 
cannot rescue the third peak, at least for |n run | < 0.2. 

We note that our results are somewhat prior depen- 
dent: the f2 c dm parameter is a derived parameter in 
the standard parametrisation used by the Cosmo-MC 
package and therefore the implied prior on it is cer- 
tainly not flat. However, we feel that priors employed 
are actually a sound set of physical priors and there- 
fore our results should be fairly insensitive to any sen- 
sible reparametrisations. 

Finally, we note that there is still a possibility that 
a version of MOND theory with a high third peak is 
constructed. However, this would seem rather arti- 
ficial. The ratio of height of third to the first peak 
encodes the amount of the baryonic drag and in or- 
der to get a high third peak one needs some cold dark 
matter like element. Even if this eventually turns out 
to be a scalar field in a MONDian theory, the present 
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data indicate that dark matter theory is at least a very 
good approximation to the full underlying theory. 
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